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ABSTRACT 


The purpose of this thesis is to analyze the vertical plane response of surface 
ships in close proximity towing. The problem is formulated by using the heave and pitch 
equations of motion in regular waves. The vertical motion of the leading and trailing ship 
attachment points is calculated. The relative motion between these points is then matched 
through a notional spring/damper model of the connection. This allows calculation of the 
complete response amplitude operators for the two ships in terms of their relative motion 
and connection force. Parametric studies are conducted in terms of connection spring and 
damper characteristics, speed, and sea direction. Regular wave results are extended in 
standard fully developed random seas. A notional example provides insight into future 


studies necessary to validate the close-proximity towing concept. 
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I INTRODUCTION 


A. PROBLEM STATEMENT 


Towing of large payloads with small, power dense vessels is a proven means to 
cost effectively transport cargo. Divorcing the prime mover from the load bearing vessel 
results in a rotich larger payload fraction dedicated to cargo. This separation also allows 
the cargo vessel to be customized for special circumstances without altering the 
configuration of towing vessel. Another advantage is the reduced cost of utilizing a 
single sensor suite in the towing vessel to provide navigation and control of various cargo 
platforms. However, use of towing vessels has been limited to low speed, high payload 
barges. Several factors contribute to the traditional prejudices against using towing 
vessels for high-speed, medium payload operations. Conventional ocean going tow rigs 
employ long lines to diminish interaction forces between tug and tow. This results in 
poor maneuverability in constrained waterways. Large interaction forces result from the 
difference in response to a seaway between the tug and barge. Thus, high magnitude 
forces due to random seas result in peak amplitudes that render towing operations 


dangerous to personnel and equipment. 


Until recently, the risks of high-speed towing have traditionally outweighed the 
rewards, leading to little interest in its development. Introduction of SWATH and related 
hull types such as SLICE that minimize sea surface interaction effects on vessels sparks 
renewed interest in the feasibility of high-speed towing based on the aforementioned 


advantages. 





B. RESEARCH APPROACH 


1. Table of Offsets Generation 


With SLICE and KAIMALINO identified as suitable platforms for study, 
background data is generated for these vessels. SLICE lines drawings and the resultant - 
table of offsets used as input into the modeling software were generated in “Seakeeping | 
Characteristics of Slice Hulls...” by Lesh in six degrees of freedom. This study verifies 
the published displacements and operating characteristics of SLICE in six degrees of 
freedom. Existing hull lines and operating environment, i.e. salt water, regular wave 
response, and motion prediction in six degrees of freedom, are verified prior to 
development of KAIMALINO data files. Next, lines drawings of KAIMALINO are 
converted by hand to a table of offsets used in response prediction. KAIMALINO 
architecture is itadiiad and verified against published characteristics provided by 
Lockheed-Martin Marietta Corporation. Table 1 summarizes this comparison and shows 
adequate agreement between published and calculated characteristics for parametric 


study of vessel response. 








Vessel Computed Published Length 
= | Displacement | Displacement | (LOA). 
KAIMALINO 


Table 1. Modeled vs. Published Characteristics 











2. Design Process 


Basic fundamentals of naval architecture are employed to determine the feasibility 
of close-proximity towing operations. A commercial FORTRAN based code, SHIPMO 
is used to model motions of SLICE and KAIMALINO in a seaway. MATLAB based 
codes are used to verify individual ship motions, calculate vessel interactions, and predict 
regular and random wave response of the integrated towing unit. Graphics, parametric | 


studies and product data are generated in the MATLAB environment. 
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Hi. SHIP MOTION MODELING 


A. OVERVIEW 





Motion of a rigid body in 3-D space can be described in 6 degrees of freedom. 


Three translational (surge, sway, heave), and three rotational (roll, pitch, yaw) 


displacements are required. The diagram below illustrates these motions. 


_z-rotation = yaw 
z-translation = heave 






) 


y-rotation = pitch 
y-translation = sway 


Vy 


x-rotation = roll 
x-translation = surge 


Figure 1. Rigid Body Model. 


Rigid body motion includes velocity and acceleration components. Throughout 


this work and in the associated computer simulations, as well as in tables and figures, the 


following symbols will be used to describe ship’s motion in a body reference frame. 





Displacement Velocity Acceleration 


1), — surge 7}, — surge, vel 1}, — surge, accel 
I, — Sway 1. — sway, vel 1], — sway, accel 
nN, — heave 1), — heave, vel 1}, — heave, accel 
1, — roll 1), — roll, vel 1, — roll, accel 
1; — pitch 1], — pitch, vel 1], — pitch, accel 
eg — yaw 1], — yaw, vel Ti, — yaw, accel 


Since both SLICE and KAIMALINO are modeled, a second subscript is added to 
each term to indicate the appropriate vessel. For SLICE, the second subscript is ‘s’ and 


for KAIMALINO the subscript is ‘k’. For example: 


1),, 7 surge, SLICE 
T,, — sway velocity, KAIMALINO 


B. MODELING A SHIP IN WAVES 


1. Background 

Ship response to head seas is “a complicated phenomena involving interactions 
between vessel] dynamics and several distinct hydrodynamic forces” (Cummins). In other 
words, an arbitrary hull form will respond to a random sea state in a random, non-linear 
manner. However, a seaway’s characteristics can be modeled and the ship’s response 
approximated as linear. In a real seaway “six non-linear, differential equations of motion 
must be set up and solved simultaneously for six unknowns” (Cummins). Advanced 
calculus and hydrodynamics theories developed by Ogilivie, Cummins, and Wehausen 
have shown that response can be reduced into a Newtonian Spring-mass-damper form that 
is frequency dependent. Further, in the case of slender hulls and moderate sea states the 


six non-linear equations reduce to two sets of three uncoupled equations. The 


6 





longitudinal motions (surge, heave, pitch) are decoupled from the transverse motions 


(sway, roll, yaw). 


Zs Frequency Dependent Equations of Motion 


A ship’s interaction with a given seaway may be modeled as a spring-mass 


damper system, in its simplest form: 


[M }i+[B]i+(ClHi=[£,] 


[M] = Mass of vessel and moments of inertia. (6x6) 
[B] = Hydrostatic damping, due to energy dissipated in wave making. (6x6) 
[C] = Restoring force and moment constants due to buoyancy. (6x6) 


[Fex] = Excitation forces and moments from seaway. 


These equations provide a basis for understanding the model; the actual equations 
of motion are slightly more complex. First, the real [M] matrix includes inertia and cross 
coupling terms [m], and is summed with an “added mass” matrix [A]. An added mass 
coefficient corresponds to each of the degrees of freedom. The added mass terms 
represent additional mass and mass moments of inertia of seawater moved when a ship 
moves in any of the six degrees of freedom. Quite simply, it accounts for the mass that 
must be supplanted by the hull as it moves in any direction. The [Fex] forces consist of 
the Froude-Kryloff and diffraction exciting forces and moments. The Froude-Kryloff 
force is due to the direct interaction between the body and the wave front. This 
interaction changes the hydrodynamic velocity profile from the undisturbed case, and 
creates the second element of the [F.,] matrix, the diffraction force. The elements of the 
[A], [B],[C], [Fex], and [m] matrices must be found to solve for the on Nk: 


f 











While the elements of the [C], [m], and [f,] (Froude-Kryloff) can be found using 
analytic equations in “Principles of Naval Architecture, Vol III,” the elements of the [A], 
[B], and [fais] (diffraction) matrices are found using “strip theory”. Strip theory involves 
examination of a two-dimensional “strip” in the x-y plane of the vessel. “The flow field 
along this strip is approximated by the assumed two-dimensional flow along the strip” 
(Cummins). To obtain the effect on the entire vessel, all the strips are integrated along the 


ship’s length. This process is repeated in each element to determine the second group of 
L 
matrix elements above. For example: A,, = { a, ax. Further discussion of strip theory and 


the elements of the [A], [B], and [C] matrices are available in “Principles of Naval 


Architecture, Vol. Hf.” 


The other matrices and a more complete description of the equations of motion 


[M]n+[8)n+[ Cc] =[F.,] 
[M] = [m+A] (6x6) 
[B] = Hydrostatic damping, due to energy dissipated in wave making. 


[C] = Restoring force and moment constants due to buoyancy. 
[Fex] = [fi + faite] (6x1) 


3. Equations of motion in the Frequency Domain 


The equations of motion discussed above in the time domain are valid for zero 
forward motion in head seas. The matrix elements are valid for any given frequency of 
waves. To more accurately predict motions in waves for a given forward speed and wave 
angle, the frequency of encounter, @, dictates the values in the matrix elements. Since 


linear theory requires that vessel response be directly proportional to wave amplitude at 
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the perceived frequency of incident waves, for regular waves, the vessel motions will be 
sinusoidal and have the form: 7, (t)=7),e"’, k=1...6 and 7,= complex amplitude of 
vessel response in the k™ direction. The frequency of waves @ must therefore be shifted 
to account for vessel speed (V), and relative direction of encounter, (B). “Regular” or 


sinusoidal waves may be described by: wave = 1, cos(kx ~ @t) 


No -- wave amplitude 
@ = 2n/T | 


k =2n/X -- wave number 


In deep water where depth is greater than (A/2), the dispersion relation states 
@=./kg . Thus for a given wavelength, wave frequency is known and the frequency of 
encounter is @. = (@ - kVcos(B)). 


Note that the time domain ODE’s are linear, and the output motions are in 


complex form. It is evident then that the ODE’s are easily transformed into the frequency 


domain, where the 77, ’s can be solved using algebraic methods. 


y=Ne™* 
7] =j on eit 
n= —@, Tye 


Since the exponential exists in all terms, it is canceled and the equations of 
motion in the frequency domain become: 


—[m+ A]ijo? +[Blio,7 +(Cii =[F.,,] 
if + A=-@?[m+ A]+[Blio, +[C] 
then > Afj =[F,,] 
and — fj =inv(A)[F,,] 
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C. SIMPLIFICATION OF EQUATIONS OF MOTION 


i. Decoupling Transverse and Longitudinal Motion 


The motions due to regular waves of a given wavelength and direction are now 
determined for a vessel with forward speed (V). Transverse and longitudinal motions are 
actually decoupled and may be solved as two distinct 3x3 systems vice the 6x6 system 
shown above. With this in mind, surge, heave and pitch responses and the resultant 
interactions between SLICE and KAIMALINO in these three degrees of freedom will be 
analyzed from this point forward. Surge motion may also be neglected because in long, 
slender ships, surge effects are small relative to heave and pitch. To simplify the 
equations of motion, all motions except n3 and 1s are set to zero. The expanded 


equations of motion in two degrees of freedom become: 


Analysis of the other degrees of freedom is possible with slight modifications of 


the matrix elements and variables of interest. 


10 





2. Data Generation and Computer Simulation 


Elements of added mass, hydrostatic damping, and force matrices in the equations 
of motion are readily calculated using a motion analysis program that employs strip 
theory. The Fortran based code “Shipmo” is used to generate the data analyzed herein. 
“Shipmo” takes a table of offsets as input and calculates motions and added mass 
coefficients for a range of speeds and wave angles. The KAIMALINO table of offsets iS 
hand generated from analysis of detailed scale drawings produced by Lockheed-Martin 
Marietta Corporation. The table of offsets is accepted along with a host of functional 
inputs such as wavelength, wave angle, forward speed, wave type, surge and roll 
damping, and the location of the position on the ship where the motion is to be analyzed. 

Ship motions and the matrix coefficients are calculated for multiple speeds and 
wave angles, for wavelengths from twenty to one thousand feet. A database of 
“Shipmo” output files is created for both SLICE and KAIMALINO for speeds from zero 
to twenty knots in one-knot increments, and for wave angles from zero (following seas) 
to one hundred eighty degrees (head seas) in five-degree increments. The utilization of 
this database to predict the connection forces on a close-proximity tow is further 


discussed in the next section. 


D. SLICE AND KAIMALINO MODELING AND INTERACTIONS 


1. Ship Motion in Pitch, Heave at Connection Point 


The integrated tow connection point on SLICE and KAIMALINO is chosen along 
the centerline, at deck height, at the aft-most point on SLICE and the forward-most point 
on KAIMALINO. Since heave and pitch are decoupled from sway and yaw, only the 
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elements associated with heave and pitch and their cross coupling elements are needed to 


solve for 13 and Ns. 


Using the Apar factoring discussed above, the equations of motion for each vessel, 


with the bar removed for simplicity reduce to: 


As3 53,5 + Ags 515,, = 13, +F, 
As3 13,5 + Ass 55,5 = Fs, — FX, 
Ass 134+ Ags 1s, = Fy, + fi 
Ass M34 + Ags Ms = E54 — FX 


f, - connection force on KAIMALINO 
f; -connection force on SLICE | 


Xsub - distance from CG to connection point 


Since f is a reaction force, f = f, = -f,. 
‘The 1’s above cannot be solved directly. However, making the following 


substitutions the motion due to the excitation force (U,; ) and the motion due to the 


connection force (Vnj) may be solved individually assuming a unit connection force. 


13,5 = Ms, +V3,,F 
Ns,5 = Us, +Vs,.F 
N34 = Usp +V3yF 
Ns4 =Hsy tVsaF 


Cramer’s rule is used to solve for ta; and Vj; in the following equations: 
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Ass Ms, + Ags sHs,. = F5,, Ag3 V3.5 + Ags Vs,, = 1 
Ass Hs, + Ass «ls, = Fs, As; V3, + Ass V5.5 =X, 
Ag; bsp + Ags pHs.x = Fx Ag3 V3.4 + Ags iV5,, = 1 
Ass ,bs. + Ass Hs. = Fy Ass V3.4 + Ass V5.4 = ~Xy 
Kj Ags j I Ass 
- Fs, Ass; ea —x; Ass, ; 
ms) Ags; Ags, j _ As; Ags; 
Ass; Ass, j As3,; Ass, j 
Ag; Fy; Ags | I 
_ Ass; Fs; — As; ; —X; 
ma Ags; Ags, j >t Ags i Abs, j 
As; , Ass, j Azz; Ass, j 














With Ltn; and vn; thus solved, the heave (13) and the pitch (js) may now be | 


determined for an arbitrary connection force. 


E. MODELING THE CONNECTION FORCE, F 


1. Absolute Motion at the Connection Point 

Having solved for the motions due to excitation ee (Unj) and the heave and 
pitch due to a unit connection force (V,j), superposition is now employed to find the 
overall heave and pitch of the two ships with a motion dependent connection force. Since 
rotational motion (pitch) will be unrestrained at the connection points, the connection 
force will be due to the difference in absolute translation of the two vessels at the 


connection point. The absolute motions at the connection points are described by: 
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E =, —N;,x, — Slice motion 


€. =N34—1s,.%, —> Kaimalino motion 


Notice the equations above cannot be solved until a connection force is supplied. 
However, because connection force is dependent on the difference in absolute motion of 
the two vessels, a theoretical relationship between the connection force and difference in | 
absolute motion must be assumed. A generic spring-damper interface is inserted and the ~ 
matching condition ee 


f=ké, —€, + c(. ~€) — time domain 
f =(k+ic\(€,-€,)> | frequency domain 


Recall that the absolute motions are functions of the heave and pitch amplitudes: 


N35 = Hs. tV3 F 
s,s = Hs. +V5 Ff 
N34 = May tVanF | 
Ns = Us tVsuF 


The following sequence describes the mathematical steps of simplifying the matching 


condition using the heave and pitch amplitudes and simplifying variables: 


if : 

(€,-€,)=a—bf — factored heave — pitch amplitudes 
and 

(, —S.) = M35 — May — M55 + HN 

then 


a= U,, — Us, —X, Us, TX, Ms, 
D=V3, —V3z —XV5, +XV5y 
next let > K =(k +ic), 
so— f =K(a-Df) 
Ka 
1+ Kb 





connection force > f = 
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The connection force is thus solved. Furthermore, the individual motions of each 


vessel and the corresponding connection force are dependent on vessel speed, seaway 


characteristics, and the spring-damper constants [7, , = fun(A, B,V,k,c)]. 


Z: Computer Modeling 
Modeling of SLICE and KAIMALINO’s response to regular waves is 


accomplished in all six degrees of freedom for a given set of input conditions as 
previously described. The purpose of this work is to further research the individual ship 
responses and develop a model that accurately predicts and if possible optimizes the 


connection force on a “hitch” connecting SLICE and KAIMALINO. 


Vessel response and matrix coefficients of the motion variables are found in 
“Shipmo”, which produces output files containing response and matrix data for reading in 
the MATLAB environment. These files are downloaded into the program 
“Samplemain” developed by Papoulias and improved by Nash to accommodate vertical 
and horizontal motions as cpus “Samplemain” repeats some of the functions of 
“Shipmo”, calculating heave and pitch response to a seaway as described in the preceding 
analytic discussions. The response is compared to the output of “Shipmo” prior to 


proceeding with matching condition and force calculations. 


With the SLICE and KAIMALINO heave and pitch response verified, matching 
condition and connection force calculations are added to determine both the connection 
force for optimization purposes, and the front and rear ship heave and pitch response with 
a rigid connection attaching the two. The program is currently capable of determining 


the coupled heave/pitch response and connection force for user supplied spring-damper 
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constants. With the existing strip theory database produced from “Shipmo” runs, this is 
possible for ranges of speed from (0-20) knots-every knot, wave angles (0-180) degrees- 


every five degrees, and is good for wavelengths from 20 feet to 1000 feet. 


F. | REGULAR WAVE RESULTS 

Formulation of the equations of absolute motion at the connection point (& and 
€.) and computer modeling discussed previously form the pillar for this research. The 
most interesting problem facing designers of a close-proximity towing system is 
engineering the connecting apparatus. The design must exhibit adequate strength to 
withstand the forces imposed by the differential motion of the two vessels. Using the 
mathematical model discussed above, and the MATLAB code "Samplemain.m", the 
absolute motion at the connection point and the resulting connection force is evaluated. 
The force is normalized for a one-foot wave height, and plotted with absolute motion 
versus wave frequency. Typical regular wave results are plotted for fifteen-knot forward 
speed and 180° wave angle (head seas). Absolute motion magnitude (xi,;= E. | xine E,. b 
and phase angle, and the connection force are plotted for different combinations of spring 
constant and damping coefficient values. With c=0, a small, medium, and large spring 
constant relative to displacement is chosen for comparison. A similar combination in 


damping coefficient is compared for k=0. 
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K=0 lbf/ft; C=0 Ibf-s/ft 
Absolute motion magnitude vs frequency (rad/sec) for k= 0 c=0 
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Figure 2. Absolute motion mag. K=0, C=0 


Absolute motion phase angle vs frequency (rad/sec) for k= 0 c=0 
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Figure 3. Absolute motion phase K=0, C=0 
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Connection force magnitude vs frequency (rad/sec) for k= 0 c=0 
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Figure 4. Force magnitude, K=0, C=0 
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K=5000 Ibf/ft; C=0 Ibf-s/ft 
Absolute motion magnitude vs frequency (rad/sec) for k= 5000 c=0 
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Figure 5. Absolute motion mag. K=5000, C=0 


Absolute motion phase angle vs frequency (rad/sec) for k= 5000 c=0 
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Figure 6. Absolute motion phase K=5000, C=0 
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x 10° Connection force magnitude vs frequency (rad/sec) for k=5000 c=0 
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Figure 7. Force magnitude, K=5000, C=0 
Connection force phase angle vs frequency (rad/sec) for k=5000 c=0 
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Figure 8. Force phase angle, K=5000, C=0 
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K=500000 Ibf/ft; C=0 Ibf-s/ft 
Absolute motion magnitude vs frequency (rad/sec) for k= 500000 c=0 
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Figure 9. Absolute motion mag. K=5x 10°, C=0 
Absolute motion phase angle vs frequency (rad/sec) for k= 500000 c=0 
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Figure 10. Absolute motion phase K=5x10°, C=0 
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Connection force magnitude vs frequency (rad/sec) for k=500000 c=0 
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Figure 11. Force magnitude K=5x10°, C=0 


Connection force phase angle vs frequency (rad/sec) for k=500000 c=0 
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Figure 12. Force phase angle K=5x10°, C=0 
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K=0 Ibf/ft; C=5000 Ibf-s/ft 


Absolute motion magnitude vs frequency (rad/sec) for k= 0 c=5000 
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Figure 13. Absolute motion mag. K=0, C=5000 


Absolute motion phase angle vs frequency (rad/sec) for k= 0 c=5000 
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Figure 14. Absolute motion phase K=0, C=5000 
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Connection force magnitude 


Connection force phase angle 





Connection force magnitude vs frequency (rad/sec) for k=O c=5000 
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Figure 15. Force magnitude K=0, C=5000 


Connection force phase angle vs frequency (rad/sec) for k=O c=5000 
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Figure 16. Force phase angle K=0, C=5000 
24 


5. 


PHASE ANGLE(i,,xi,) 











K=0 Ibf/ft; C=500000 Ibf-s/ft 
Absolute motion magnitude vs frequency (rad/sec) for k= 0 c=500000 
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Figure 17. Absolute motion mag. K=0, C=5x10° 


Absolute motion phase angle vs frequency (rad/sec) for k= 0 c=500000 
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Figure 18. Absolute motion phase K=0, C=5x10° 
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Connection force magnitude vs frequency (rad/sec) for k=O c=500000 


12000 


10000 


8000 


f°?) 
oO 
(2) 
Oo 


Connection force magnitude 


ih 
© 
© 
few ] 


2000 





omega 


Figure 19. Force magnitude K=0, C=5x10° 


Connection force phase angle vs frequency (rad/sec) for k=O c=500000 
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Figure 20. Force phase angle K=0, C=5x10° 
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=5000 Ibf/ft; C=5000 Ibf-s/ft 
Absolute motion magnitude vs frequency (rad/sec) for k= 5000 c=5000 
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Figure 21. Absolute motion mag. K=5000, C=5000 


Absolute motion phase angle vs frequency (rad/sec) for k= 5000 c=5000 


PHASE ANGLE(xi,,xi,) 





0 0.5 1 1.5 2 2.5 3 3.5 
omega 


Figure 22. Absolute motion phase K=5000, C=5000 


2 


12000 


10000 


Connection force magnitude 


Connection force phase angle 





Connection force magnitude vs frequency (rad/sec) for k=5000 c=5000 
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Figure 23. Force magnitude K=5000, C=5000 


Connection force phase angle vs frequency (rad/sec) for k=5000 c=5000 


50 


© 


on 
o 


ak 
=) 
o 


=t 
o 
-) 





-200 
0 0.5 1 1.5 2 2.5 3 3.5 


omega 


Figure 24. Force phase angle K=5000, C=5000 
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7. Regular Wave Results Observations 


Several fundamental properties of the system are evident in the preceding plots. 
Figures 2 and 3 show the vessel response for (k=0, c=0). This response is equivalent to 
two independent vessels in single file, i.e. disconnected. The offset in @ between &, and 
€ shown in figure (2) is due to the distance separating the connection points. The 
absolute adios shows resonant peaks for A = 4-6 times ship length at the example speed 
and heading: Figure (13) shows the magnitude of these peaks is cut in half as damping is 
added (c=5000 Ib-s/ft). Such peaks evident in regular waves should be substantially 
reduced in a real seaway because of the random nature of real waves. Furthermore, as 
would be expected in the disconnected case, figure (4) shows zero connection force for 


all @ when the spring constant and damping constant are zero. 


For average displacement of SLICE and KAIMALINO (~200 tons), as the spring 
constant is raised to large values relative to the displacement (k=500,000 Ib/ft ~ 223 
ton/ft), the response of the two vessels approaches that of a suale rigid vessel. This 
behavior is evident in figures (9) and (10), where magnitude and phase angle of the two 
vessels merge to the same values over the range of @. A very large damping constant 
value (c=500,000 Ib-s/ft) also models an inflexible connection as shown in figures (17) 


and (18). 


Force and absolute motion magnitude and phase angle, are evaluated to verify the 
software provides reasonable results. With results verified, the response to any given sea 
way is now known. Mapping vessel response to seaway motions is the ultimate goal of 
regular wave models. The function that maps wave input to ship response is the 


"Response Amplitude Operator", and will be further addressed in random wave analysis. 
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Plots also provide insight into the variation of force and absolute motion as spring 
and damping constants are varied. Regular wave modeling program and results are 
useful tools in the design spiral of an actual towing mechanism. Design trade-offs 
between minimizing relative motion (large k and c values), and minimizing connection 
force/tow bar size can be roughly evaluated. Regular wave modeling results are not 
however, precise enough to base actual design upon. The regular nature of the sinusoidal | 
sea can show false resonance, and abnormally high peak magnitudes that would not be 
encountered in random seas. To more precisely model vessel response, the response 


amplitude operator must be mapped to an applicable random sea spectrum. 


30 











lil, RANDOM WAVE RESPONSE 


A. BACKGROUND 


1. Spectrum Selection Criteria 


A seaway's "spectrum" is a probabilistic function developed by taking the Fourier 
transform of the correlation function for free surface elevation. (Cummins) The 
correlation function contains wave height and period data from sources such as buoy 
observations. The spectrum - S(q@) — is a measure of the energy contained within a wave 


system. In a plot of S(@) vs. @, the area under the curve represents the mean energy 


stored in a particular wave system, E= [s (@)do. 
0 


Numerous wave spectra are available as input for modeling a vessel's operating 
response to a seaway. Selecting the appropriate spectrum should be accomplished with 
due regard for the ship's expected operating environment. Environmental conditions such 
as wind and swell vary geographically, and ship design should be tailored so the vessel 
responds optimally to the prevailing conditions. Since the environment that SLICE and 
KAIMALINO will operate is unknown, the Pierson-Moskowitz spectrum is chosen. This 
model predicts the wave spectrum for fully developed, long-crested seas with no 
underlying swell. Fully developed seas contain waves at equilibrium, independent of 
fetch and duration of wind. Long crested seas have parallel crests and are assumed to be 


unidirectional. 
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The Pierson-Moskowitz spectrum is described by 


@ Uw 


where: g = gravitational constant 
@ = wave frequency (rad/sec) 


U = Wind speed at19.5 m above free surface 


The above spectrum is dependent on wave frequency and wind speed, a metric the 
underlying regular wave research does not provide or account for. Regular wave results 
provide spectral vessel response as a function of frequency for given significant wave 
height. Correlation between significant wave hiotit and wind speed has been extensively 
developed, and frequency dependent spectral formulations derived based on significant 
wave height. (McCreight) Using an empirical relationship between wind speed and wave 


height, the Pierson-Moskowiiz spectrum can be predicted using the following 





relationship: 
-3 
S(@) = eels ‘ 
@ 
where: g = gravitational const. (32.2 ft/s*2) 
@ = wave frequency (rad/sec) 


Hy3 = Significant wave height (ft.) 


Significant wave height is defined as the average of the highest one-third of all wave 


height observations. 
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2. Response Amplitude Operator (RAO) 


Also known as the motion transfer function, the RAO maps the complex response 


of a vessel to a seaway or input spectrum as a function of frequency. 
S,(@) =|RAO(@)|” S(@) Where Sp(@) is the response of the vessel to the input sea 


spectrum for a given frequency. This very powerful relationship allows motions and 
terms derived from ship motions to be predicted for a given wave frequency and 
significant wave height. For example, the complex absolute motions predicted in regular 


wave modeling (€,, &,) are converted into RAO’s for absolution motion: 
RAOE,,) = abs(E,,)- 
The response spectrum for absolute motion is: 
Sp —€,,(@) =|abs(E,,)[ S(@). 


This process may be applied to all motions, and in the case of the close-proximity towing 


system, the connection force response is: 


* S(@). 





Sr _— dT senstien (@) = labs( f connection ) 


With the spectral response of a vessel’s motion thus determined, the design spiral 
continues, with random wave results providing a more complete assessment of design 
objectives. In order to conduct trade-off analysis, or to evaluate performance against 
changes in environmental, operational, and design parameters, the statistical properties of 
the response must be determined. In other words, while it is useful to predict the 
response for a given eae frequency, it is more productive to compare performance over 
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the entire range of the spectrum. A good measure of performance over a range of 
frequencies is the significant double amplitude of the response. Double amplitudes are 
obtained by simply integrating the response with respect to frequency over the frequency 


range of the input spectrum. For instance, the significant double amplitude of the 


@ 
ig 
absolute motion of SLICE at the connection point is: O, = | S,_6,(@)do. 


Wo 


B. RANDOM WAVE MODELING OF SLICE AND KAIMALINO 


1. Process 


Regular wave modeling of the SLICE-KAIMALINO integrated tow rig discussed 
in the previous chapter yields RAO's for absolute motion of both vessels, as well as the 
RAO for connection force. Using variable forward speeds and wave angles yields RAO's 
that are functions of the frequency of encounter @, rather than actual wave frequency @. 
S(@) is readily converted to S(@) because the energy of the seaway will en constant 


whether viewed from a stationary point or a moving ship. 





S(@)d@=S(@,)d@, .. {energy(@) =energy(Q, )} 
da a’ 

S(@) = S(@,)— {o@, =@-—U cos B} 
da@ g 


S(@,) =S(@)[1 7 Ucosp]" 


The first operation performed in the random wave analysis software is defining 
the Pierson-Moskowitz spectrum and transforming it as shown above. Next, the response 


spectra are defined. The most interesting response in determination of feasibility of the 
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close-proximity towing system is the connection force, whose response spectrum is 





2 
defined as S,(O,) =|abS( fronection)| S(®,)- 
Ogg 
Integration of the connection force response spectrum, O gonn = | 5 ,(Q, )da, 
Deo 


is accomplished numerically by summing the product in the integrand of the preceding 


integral. In other words, 


#of@'s 
O =Srot DY Spat Sy @. gy Pe) 


fconn 


i= 


The resultant Ofonn 18 now transformed to significant double amplitude, 
o, =4, lo fon - Lhe significant double amplitude represents the average of one third of 


the highest probable connection forces encountered for the given input condition in 


waves with wavelength from twenty to one thousand feet. 


2: Results 


The random wave simulation described above adds significant wave height to the 
list of input parameters that were varied in the regular wave studies. Recall that a 
database of regular wave RAO's for both SLICE and KAIMALINO was created, for 
speeds from zero to 20 knots, and wave angles from 0° to 180°. Combining the database, 
regular wave simulation, and random wave simulation enables parametric studies to be 
conducted. Several questions should be answered before the tow mechanism is designed. 
Is the connection force in seas up to sea state five small enough to make integrated 
connection feasible? If the force is manageable, what spring and damping constants 


ee) 








should be used to minimize the force? And finally, which sea directions and ship speeds 


drive connection force to unacceptably large values? 


Armed with random seas software, these questions are researched by varying ship 
speed, wave angle, significant wave height, and spring-damper constants. Standard 
values are used throughout to allow cross-reference between plots. Standard speed is 15 
knots, wave angle is 45°, significant wave height is 5 feet, and spring and damper 
constants are set to zero. Similarly, when each parameter is varied, it must be done so in 


a like manner from one run to the next. Standard parameter variations are: 


1—100,000 Ib-s/ft 
K in lb/ft 
Table 2. Parametric Variations 





Using these parameter variations, random wave simulations are run and 
connection force is calculated as a significant double amplitude (6). Connection force is 
selected as the common metric against which all variables are compared. Similar 
parametric studies can be run using icles motion (€;, &) as the dependent comparison 
variable. Plotting Of versus V, B, k, and c reveals optimum values of each of the 
variables for minimizing significant connection force. Finally, the plots provide a tool 
for design of the actual close-proximity towing mechanism. The optimum spring-damper 
values and maximum expected connection force output by the random wave simulation 
provide the basis for solid mechanics engineering of the tow bar. For instance, with 
maximum 6; as the design force, maximum yield stress and Euler buckling theories might 
dictate the cross-sectional parameters of the tow bar. Such a case, study is presented in 


the next chapter. 
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a. O7 (op versus Hy1/3)-f 
x10° Sig. force vs. Sig. wave height for k= 5000 c=5000 V= 15 kts 
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Figure 25. Force vs. Hi3), (B varied) 


x10° Sig. force vs. Sig. wave height for k= 5000 c=5000B= 45° 
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Figure 26. Force vs. H 3); (V varied) 
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Figure 27. Force vs. Hi13), (kK varied) . 


5 x 10° Sig. force vs. Sig. wave height for k=0 V= 15 kts B=45° 
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Figure 28. Force vs. H(1/3), (c varied) 
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b. OF (bf) VETSUS Vues 


Sig. force vs. Velocity for H..= 5 ft k=5000 c=5000 


x 10° 
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Figure 29. Force vs. Vis, (B varied) 





Figure 30. Force vs. Vitis, (Hi3 varied) 
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Figure 31. Force vs. Vits, (k varied) 


4 Sig. force vs. Velocity for B= 45° H, ,= 5 ft k=O 
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Figure 32. Force vs. Vis, (c varied) 
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Sig. force vs. B for H,,.= 5ft V=15kts 
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Figure 35. Force vs. Baegrees, (kK varied) 


Sig. force vs. B for H, ,= Sft V=15kts k=0 
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Figure 36. Force vs. Baegrees, (C varied) 


42 







d. OF ip versus Kingp: 


4 0° Sig. force vs. Spring Const. for c=0 V= 15kts H, = 5 tt. 
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Figure 37. Force vs. kip, (B varied) 
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Figure 38. Force vs. kj, (V varied) 
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Figure 39. Force vs. Kye, (H(1/3) varied) 


Sig. force vs. Spring Const. for V=15kts B= 45° H, ,= 5ft 
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Figure 40. Force vs. King, (c varied) 
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OF (ap Versus Cipp-sif: 


Sig. force vs. Damping Const. for k=O V= 15 kts H, = St 
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Figure 41. Force vs. cip-sp, (B varied) 
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Figure 42. Force vs. Cip-s, (V varied) 
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x10’ Sig. force vs. Damping Const. for V=15kts B= 45° k= 0 
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Figure 43. Force vs. Cyp.sg, (H(1/3) varied) 


Sig. force vs. Damping Const. for V=15kts B= 45° H, jg= 5 tt 
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Figure 44. Force vs. cys, (k varied) 
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3. Seakeeping Evaluations 


Variation in connection fies double amplitude with significant wave height 
provides largely intuitive results. Figure 27 shows a linear rise in connection force as 
significant wave height increases from zero to ten feet. Increasing the Sonne constant 
values from | Ib/ft to 100,000 lb/ft is equivalent to raising the rigidity of the connection 
from disconnected to rigidly connected. As expected, the connection force rises as the 
connection becomes more rigid. The slope of the o¢ versus H,y/3) increases as spring 
constant values increase. A similar sia is evident in figure 28 where connection 
force amplitude rises with wave height, for all values of damping constant, and the slope 


of the rise increases with damping constant increase. 


Forward speed and wave angle variations reveal several interesting results. 
Figure 25 shows higher connection force for following seas (B<90°) than head seas for 
wave height less than four feet. However, as significant wave height increases from four 
to ten feet, the slope of the head seas cases Siete rises, yielding much higher 
connection force for the head seas case at a ten foot wave height. The most interesting 
results shown in figure 26, the of Sein Has) plot for various speeds are the high 
connection forces corresponding to speeds of three, six, and twelve knots. Six knots 
yielded the highest connection force, with twelve and three knots being next in line. The 
effects of speed and wave angle on connection force are reaffirmed in the 6; versus V and 
Of versus 8 plots. ae 30, 31, and 32 clearly show the six and twelve knot force 
peaks. Figure 34 confirms the much higher connection force for following than head 


seas. It is important to note however, that the connection force used in analysis was 
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developed from vertical plane motions only. Transforming the simulation to include 


surge forces is likely to raise the magnitude of the connection force in following seas. 


The relationship between spring-damper constants and connection force is 
perhaps the most useful in design of the towing mechanism. Force dependence on wave 
height, wae angle, and speed provide operating characteristics of the integrated vessels, 
and result in engineering and operating limits due to environmental factors. For the 
design in question, the only remaining variables that can be manipulated by engineers to 
minimize the connection force are the spring-damper constants. Evaluation of the 
parametric plots of Of versus k and oO¢ versus c reveals optimum spring and damping 
constants to minimize connection force. For instance, figures 40 and 44 show that for 
values of spring and damper constants in the two to three thousand Ib/ft and lb-s/ft range 
the connection force drops to values less than 8000 pounds, while providing adequate 
rigidity. Follow on structural design of the tow member can be accomplished by iterating 
the spring-damper constants in this region for significant connection force and absolute 


motion to meet designer specifications. | 
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IV. SAMPLE TOWING DESIGN 


A. NOTIONAL ARCHITECTURE 


Several notional towing mechanisms have been introduced by Lockheed-Martin 
Marietta for possible design in the integrated SLICE-KAIMALINO project, and several 
features are common to each design. As discussed in chapter I, forces on the tow rig are 
due to the difference in motion of the two vessels. A hybrid design was developed for the 
NPS Total Ship Systems Engineering SEA-LANCE high-speed patrol craft and grid 
deployment module. This design snenpe to minimize the number of degrees of freedom 
constrained, while also simplifying control architecture and tow mating in the open 


OCCal}n}. 


Figure 45 shows a notional close-proximity towing design. This design 
minimizes connection forces by constraining only those degrees of freedom necessary to 
provide control and stability in adverse sea states. The most severe motions in a sea way 
are expected to be in the form of roll, pitch and yaw. To minimize handling equipment 
size these motions are unconstrained between the SLICE and KAIMALINO in the tow 
bar. Yaw is constrained at the bow of the KAIMALINO only by "moment cables that 
prevent jackknifing. Surge is constrained by the tow bar, while sway is limited by the 
directional stability of KAIMALINO's SWATH hull and constant tension winches that 
could be mounted at the forward outermost edges of KAIMALINO's bow. Hinges that 
decouple pitch at both the KAIMALINO and SLICE extremities minimize heave forces. 
Finally, roll is decoupled on KAIMALINO and SLICE by a roll bearing at the stern 


of SLICE. 
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Figure 45. Notional Tow Connection 


B. DESIGN APPROACH AND RESULTS 


1. Operating Environment and Assumptions 

Design of a notional towing system is accomplished by using simulated 
connection force outputs to analytically calculate stresses at critical locations in design. 
Geometric and force magnitude considerations dictate tow bar length. Several stress 
analysis techniques will then be used to determine the minimum cross-sectional size of 


load bearing components. 
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The design approach and outcome are heavily influenced by simulation and 
analytical limitations of the research. For instance, only vertical absolute motions are 
simulated, so sizing of the tow bar and equipment is based on peak vertical forces 
expected at a transit speed of 15 knots. Peak vertical connection forces were discovered 
in the previous chapter to exist in head seas. The operating environment evaluated 


(B=180°), is expected to result in the highest connection forces for a given sea state at 15 


knots. 


Next, the spring-damper constants optimized in the random wave analysis must be 
modeled to predict the interaction between SLICE and KAIMALINO. Recall the 


equations of absolute motion describing the connection force 


FHhE—&) +6, -€,) time domain 

f =(K+icXG, -G,) — frequency domain 
The damping constants are difficult to model without icnawieass of the exact 
characteristics of the joint rotations and are assumed to be zero. The vertical connection 
force however is controlled by the tension between the vessels and the amount of vertical 
displacement separating them. As such, the spring constant value may be modeled as 
jo | eee where T is tension in the tow bar due - hull resistance obtained from resistance 
versus speed curves of KAIMALINO, and Liow is the length of the tow bar. Notice that 
the spring constant decreases with tow bar length. Geometric considerations dictate that 
tow bar must be long enough to prevent impact of SLICE and KAIMALINO during 
maneuver. While connection forces drop as length increases, bar rigidity decreases and is 
more prone to buckling. Additionally, bar length should be minimized to increase 


integrated towing maneuverability. 
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Ze Random Sea Modeling 


Random seas modeling developed in the previous chapter is the basis for 
connection force determination and subsequent handling gear sizing requirements. 
Spring constants in the modeling software are set to k=T/Liow, and connection forces are 
evaluated for lengths from 10 to 20 feet. For a forward speed of 15 knots, KAIMALINO 
resistance (T) is 35,000 pounds. Figure 46 shows the rise in connection force as 


significant wave height increases and tow bar length decreases. 


Sig. force vs. H, ,, for B= 180° V= 15kts c=0 


x 10° 





Figure 46. Of VS. Ha); (k=T/Liow) 
To minimize connection force and allow rotational freedom up to 45°, the tow bar 


length is chosen as Lipw=20 feet. As a result, for 10-foot seas, the maximum vertical 


force on the tow bar will be 79,000 pounds. 
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3. Stress Evaluation and Component Sizing 


Assumed forces include: forces from seaway and hydrodynamic resistance. Each 
of these forces results in a stress on the tow system. Three structural limitations are 


considered; Euler buckling, static yield stress, and shear yield stress. 


Seaway forces are derived from strip theory for a given tow bar length. The 
primary forces of concern are the vertical force applied to the tow bar both in 


compression and tension, and the towing resistance. F,cos(@)=T fF, is the axial 


resultant force in the tow bar, and @ is the angle formed due absolute vertical motion 
between SLICE and KAIMALINO. For a max expected pitch angle of 25°, and towing 
resistance of 35,000 pounds at 15 knots, F=38,618 pounds. For an assumed box beam 
with outer diameter of 8 inches, F, is used to calculate tow bar thickness, (t=0.34 inches) 


using basic buckling and static yield stress analysis as follows. 





2 
gaat NEI - ] 4 4 
Buckling: F( safetyfactor) _ 7? [= 19 Se =—s. ) 
Yield stress: F, ( safetyfactor) = 0, (s? = s’) 
Where: Esteel=29,000 psi Le = Liow = 20ft 
safety factor = 5 Oy =36,000 psi 


ee 


Table 3. Box Beam Requirements 
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The final consideration in the notional design of figure 45 is the design of the pins 
used at the pivot points of the mechanism. Shear stress is the primary concern at these 
points. The force acting on the pins is assumed to be the resultant of hydrodynamic 
resistance tension and maximum vertical force read from figure 46. 


Shear force v=./T’ + f2,, =87,934 pounds. From traditional solid mechanics, 


T ,=0.50, 
v(safetyfactor) =2T , A 


Resulting in a solid circular pin of diameter 3.94" = 4", 
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C. CONCLUSIONS AND RECOMENDATIONS 


1. Conclusions 


SLICE-KAIMALINO close proximity towing operations are feasible based on 
analysis of regular and random wave vertical plane vessel response. The close proximity- 
towing concept promises to be a cost effective means to transport a wide range of 
payload configurations at high speed. It may also result in the development of versatile 
warships capable of multiple roles as fighting ships and payload delivery platforms. The 
goal of this research, to provide an estimate of SWATH vessel motions and connection 
forces based on a generic connection has been accomplished in heave and pitch. Since 
heave and pitch are expected to be the most violent motions constrained between the 
vessels, these motions are likely to result in the highest magnitude connection forces. 
With this in mind, the analysis shows that connection forces are manageable with 
reasonably sized handling equipment. The research also reveals trends in the operating 
characteristics of the vessels, and insight into the optimization of spring-damper values 
that should be designed into the connection. In particular, head seas provided the highest 
magnitude forces. Three, six, and twelve knots yielded peak absolute motions and 
connection forces. Connection force response characteristics changed from following 
seas dominated at wave height less than four feet, and head seas dominated greater than 
four feet. A what if analysis of spring and damping constants yielded optimum values of 
k and c in the three to seven thousand range. Finally, a simplified solid mechanics 
evaluation of the vertical forces resulted in a 20 feet long box beam with side length of 8” 


and a thickness of 0.34”, whose thickness was dictated by yield stress. 
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2. Recommendations 


This thesis provides a solid foundation upon which future study of close 
proximity towing operations can be based. Data files describing SLICE and 
KAIMALINO operating characteristics have been developed for motion in six degrees of 
freedom over a large range of environmental conditions. Related software and the 
process i analyze random wave seenits and design a notional tow bar are outlined. 
However, several important follow on studies should be conducted to fully validate the 
Close proximity-towing concept. First, the software should be modified to include 
absolute motions and connection forces in all Six degrees of freedom to confirm the 
assumption that vertical forces will be most significant. Next, further analysis and more 
complete modeling of the spring and damper constants as they relate to a nacional: 
Structure should be conducted. A finite element model of the connection should be 
constructed to fully evaluate stress states and critical locations. Finally, tow tank scale 
models of the notional design should be built and evaluated to confirm simulation results 


and determine the effects of close field interactions between the vessels. 
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APPENDIX A 


Input file, SHIPMO.IN for running regular and irregular wave analyses on the SLICE 
hull form. Refer to Appendix A of the SHIPMO.BM User’s Manual for format and line 
content information. 


SLICE HULL FORM GENERATED BY D.B. LESH APRIL 1995 
Updated by NASH JAN 2001 
Vertical and horizontal motions with updated surge 


damping 

0 0 

1 0 1 
105.0000 
1.6557E+02 
33.0000 -26. 

1 48.6000 
16.5000 

5 44.8750 
16.0000 
16.2500 -0Q. 
16.5000 ~1. 
16.7500 -Q. 
17.0000 0 

8 40.8750 
15.5000 0 
15.8000 -2 
16.1000 -4 
16.1000 ~-10 
16.9000 -10 
16.9000 ~4 
17.2000 -2 
17.5000 0 
L5 39.8750 
15.4000 0 
15.6000 -1 
15.5500 -7 
14.6250 -8 
14.5000 -10 
14.6250 -11. 
15.0850 -11. 
16.5000 -12 
17.9000 -11 
18.3750 -11 
18.5000 -10 
18.3750 -8 


QO. 


QO. 


0 
20 


0 
1 


1.9905 
0.0000 


0000 

0 
0000 

0 
0000 
9000 
8000 
9000 


.0000 


0 


.0000 
.0000 
.0000 
.0000 
.0000 
.0000 
.0000 
.0000 


0 


.0000 
.5000 
.3100 
.9200 
.0000 


1000 
4100 


.0000 
-4100 
. 1000 
.0000 
-9200 


0 1 0 0 0 


32.1740 1.26E-05 


1.0000 


.0000 


.0000 


.0000 


0000 


59 


0 


0 





17 
He 
17 
iB 


15: 
.0000 
.5800 
.3000 
14. 
.5000 
.-2900 
. 7000 
-4200 
18. 
.0000 


a5 
13 
13 


16 
18 
Lo 
19 


18 
i> 


14. 
14. 
14. 


13 
12 
i3 


16 


20 


18 


13 
pe Bs 
15 
13 
12 


16 


18. 
20. 
.5000 
20. 
.8000 
. 8000 


20 


17 
177 
13 
15 


4500 
4000 
.6000 


37 
0000 


7000 


0000 


33 
8750 
8750 
8500 


.0400 
.5000 
.0400 
14. 
. 9000 
18. 


5000 


5000 


.0000 
20. 
20. 


5000 
0000 


-1500 
iG 
18. 


1250 
1250 
23 


. 2000 
. 2000 
.0400 
.5000 
132 
14. 


0400 
5000 


.5000 


5000 
0000 


0000 


19 


.8750 


ele 


. 0000 


.0000 


.0000 


- 0000 


60 


0 


0 


0 


0 





QO. 


.9100 -6 
6640 -6 
.0380 -10 
.5000 -11 
100 13 
.5000 -13 
.2300 -13 
.5000 -11 
.0000 -10 
.3400 -6 
.1000 -6 
.1250 0 
16.8750 
.4000 0 
4000 -7 
.0600 -7 
0100 -8 
.6200 -10 
.0100 -11 
0600 -12 
.5000 - -12 
.0000 -12 
.1000 -11 
.3800 -10 
.1000 -8 
20000: ~=7 
6000 -7 
.6000 0 
7.1250. 
5000 -9. 
1400 -9. 
8800 -9. 
.7800 -10. 
.8800 -10. 
.1400 -10. 
.5000 -10. 
.8600 -10. 
51200: =10: 
2200: =10: 
1200 -9. 
8600 -9. 
5000 -9. 
0.0000 
.5000 
-10.1250 
.5000 -10. 
-11.3000 


0.0000 


0.0000 


0 


1 








23 
22 


22 
23 


25 
25 
25 


23 
13 
23 


22. 
Zig 
20% 
Zi, 
.1500 
.5000 
.0000 
26. 
-4000 
26. 
25%% 


22 
23 
25 


26 
23 
15 


22 
23 


19 


20. 
-5000 
.9000 
.5000 
2s 
.5000 
.0000 
.5000 


21 
25 
25 


27 
27 
25 


23. 
24. 


11 
22 
22 
21 


.2000 
.4000 
21. 
21. 
213 


5900 
3000 
5900 


.4000 
.9000 
24. 
.4000 
. 7000 
.4000 
24. 
.9000 


6000 


6000 


-13 


.5000 


1000 
0800 
6000 
1800 


0000 


0000 
0000 


-5000 


-16 


.4000 
-4000 
21. 
20. 
.5000 


5000 
0000 


0000 


0000 


6000 
6000 
-19 


. 1000 
- 1000 
- 3400 


- 3500 


.0000 


.0000 


.0000 


62 


1 


0 


0 





19: 
20. 
eo HA 
26. 
21% 
22% 
24. 
24. 


13 


2A 
oe 
214 
20. 
192 
20 :« 
23° 
26: 
vd a 
27. 
Zoe 
24. 
~ 24. 


i3 


22% 
ade 
22% 
Zi 
20°. 
2A: 
23. 
ADs 
20s 
i Ao 
aos 
ao. 
Pio 


15 


a3 
2o% 
23. 
22. 
22% 
22% 
a3 
Ae Te 
24. 
24. 


= 14. 


age Oe 
=10% 


.0000 


.0000 


.0000 


63 
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e) 


0 








24.6800 -10 
24.5200 =9 
24.0000 -8 
23.6000 -8 
23.6000 0 
1 -46.1250 
23.5000 -10 
0.1000 0 
2.0400 
2962 
-46.00 0 
1.0000. 
00.0000 
0.1000 
005.0000 005. 
0.0 


.0000 
-4100 
£97350 
.8000 
.0000 


0 


.0000 
.0000 


0.0 


.0000 


.0000 1 


9.5000 


20.0000 1000.0000 
0.0000 


0000 


0.0000 


20.0000 


00.00 


APPENDIX B 


' Input file, SHIPMO.IN for running regular and irregular wave analyses on the 
KAIMALINON hull form. Refer to Appendix A of the SHIPMO.BM User’s Manual for 
format and line content information. 


KAIMALINO horizontal and vertical motions 





With surge damping 
Generated by C.A. Nash Jan 2001 


0 0 0 0 1 0 0 0 
1 ) 1 20 1 
80.50000 139905 32.1740 1.26E-05 
2.6500E+02 0.0000 

33.0000 -26.0000 1.0000 
15 40.5000 0.0000 a 
19.6120 ~-11.0630 

1932930 =11.2630 

19.1560 -11.4370 

19.0900 -12.0000 

19.1900 -12.3800 

19.3700 -12.6000 

19.7400 -13.0140 

20.1980 -12.9950 

20.5630 -12.8440 

20.8440 -12.5630 

20.9950 -12-..1980 

20.9950 -11.8020 

20.8440 -11.4370 

209930 =1L1.1560 

20.1980 -11.0630 

3 38.5000 0.0000 1 
19.4100 -9.0100 

18.3100 -~9.4600 

17.4700 -10.3100 

16.9700 -11.7000 

17.8500 -14.1500 

19.1200 -14.9200 

20.3000 -15.0300 

21.6900 -14.5300 

22.5400 -13.6900 

23.0340 -12.2990 

22.9200 -11.1100 

22.1600 -9.8400 
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20. 


15 


18. 
.5300 
.5300 
-2000 


19 
9 
18 


A ha a 
16. 
4500 « 
19. 
. 3200 
.6900 
- 2400 


17 


—20 
22 
23 


22. 


20 
20 
21 
14 
18 
19 
18 
- 16 
16 
17 
19 
ZA. 
22 
Zo 
23 
22 
20 


14 


18. 
.8800 
-2500 
. 7100 
.0800 
.1600 
.0000 


18 
17 
16 
17 
18 
20 


21. 
. 1700 
- 2900 
. 3400 
- 3100 
. L000 


23 
23 
22 
21 
20 


21. 


8900 
34 
5200 


5000 
7800 


0000 


6400 


.4700 
.4700 
-2100 


28 


.2000 
.0700 
.1800 
. 7400 
.8600 
.6800 
. 3600 
. 2600 
. 5400 
.2600 
. 1450 
. 3200 
-9400 
24% 


8100 
26 
0200 


8400 


97700 


.9000 


252 


0000 


.0000 
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0 


0 


0 





0. 
.6700 


22.5000 
.2/00 -Q. 
.0000 -8. 
.8250 —o 
.2000 -10 
.6400 -11. 
20300" = 1S 
FLO Sia. 
13300: =15 
-5900. =<i24, 
.8000 -13 
(3600: sl. 
.6100 -9 
.2800 -8. 
. 7000 =O. 

20.5000 
.4800 
.9900 -7 
.8600 aaa an 
.2300 =o 
.9000 -11. 
-2500: =14. 
-6800 -16. 
.4500 -16. 
-5200 -15 
.7700 -14. 
.5400 -12. 
.7700 =o 
-2200 -7 
.5100 =(Q'2 

15.0000 
.7700 a oa 
.0000 -7 
.9000 rc ge 
.0300 2 
i2900° -=11-. 
.4300 -13 
s9700° =13: 
#0700 -=16. 
.4700 -16. 
.6500 -15 
-6800 -12 
.6800 -11. 
.6900 -8. 
.2700 =} 
.2900 =O). 

6.5000 


0 
2500 
6300 
8500 


.1300 


6700 


.9900 


1200 


. 3560 


9700 


.8700 


6700 


.8600 


8900 
2500 

0 
0000 


8000 
4800 
1100 
5200 
3400 
5200 


st730 


5200 
0000 


-4800 
.8400 


0300 
0 
8500 


-2200 


1500 


- 3500 


5300 


. 3900 


6900 
6800 
7500 


.9700 
.9300 


0700 
9700 


-2400 


9500 


.0000 


.0000 


.0000 


0.0000 


67 


1 


0 


0 


0 








0.0000 1 


0.0000 1 


0.0000 1 








-3400 
-0200 
.2300 
.5700 
. 3400 
.0700 
.0000 
.9400 
.6600 
.4300 
. 7800 
.9900 
.6600 
.4700 


-7 


.9300 
.3400 
.0200 
.2300 
.5700 
. 3400 
.0700 
.0000 
. 9400 
.6600 
.4300 
. 7800 
.9900 
.6600 
.4700 


-11 


-5500 
.0300 
-5300 
.0200 
.2300 
.1300 
. 3400 
.0700 
.9400 
.6600 
.4300 
. 7800 
.9900 
. 2400 
. 3700 


-15 


2a: 
~9. 
-i1. 
-13. 
-15. 
-16. 
-16. 
-~16. 
-15. 
213: 
~11. 
-9. 
8, 
27 


.5000 


-7. 
sa = 
_9. 
a nee 
= 
=: 
= 16% 
= 16:2 
=O. 
=o 
= 13% 
~11. 
~9. 
_3. 
oe 


.5000 


-~0O. 
-1. 
-7. 
=o 

=<11. 

-14. 

“15. 

-16. 

-16. 

Ses 

=1 3) 

-i1. 
aos 
cal ae 
-Q. 


5000 


0.0000 


0.0000 





dl 


0 


0 








18. 
.1700 
.1900 
-4800 
.8100 
-4600 
-1200 
.0000 
- 1400 
24. 
24. 
.8800 
.8100 
.8300 
.2700 


19 
19 
1'7 
15 
15 
7 
20 
22 


22 
20 
20 
21 
15 
18 
18 


18 
16 
16 
i 
20 
22 
23 
23 
23 
zd 
21 


15 
18 
18 


18 
iB 


1 i 
.0000 
ZL. 


20 


22 
22 
22 
pa 
21 


7200 


5100 
1900 


= 19 


.2400 
.4500 
20. 
.6200 
.6400 
.2100 
. 3000 
.0000 
.0000 
. 1700 
.8000 
- 3600 
. 3800 
-59600 
Zs 


0000 


7600 
~23 


.1300 
.6200 
18. 
.6200 
.0100 


1400 


8800 


1500 


.59000 
9900 
. 3200 
. 3800 
. 3800 
Zi. 
dics 


8600 
8700 


0. 
=2 
-7. 
-8. 

-10. 
=i 
=15 
-16. 
ee oe 
seat Oe 
=1.0) 
=O:5 
-l. 
a2. 
-Q. 


.5000 


=O... 
=o 
-8. 
-8. 
=10:, 
=e. 
-14. 
= is 
aa Bs 
-14. 
gl 
-10. 
—O:. 
Hes 
-Q. 


.5000 


=. 
= 
= é 
29. 
a ee 
-14. 
+153 
-14. 
=13. 
= Ll 
-10. 
2O.. 
=>. 
oD. 
-Q. 


.0000 


.0000 
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0 


14 -32.7500 0.0000 0 
19.1800 -0.6500 
18.7600 -2.9700 
18.7600 -5.9700 
20.0000 -11.2700 
19.4400 -11.5400 
19.2800 -12.0000 
19.5400 -12.5600 
19.8600 -12.7100 
20.4000 -12.6000 
20.7300 -12.0000 
20.4600 -11.4400 
21.2500 -5.9700 
21.2500 -2.9700 
20.8200 -0.6500 
11 -36.0000 . 0.0000 1 
19.8900 -2.2500 
19.6900 -5.9700 
19.5000 -6.1500 
19.6100 -9.9100 
19.8200 -13.7900 
19.9700 -14.6200 
20.1800 -13.7900 
20.4000 -9.9100 
20.5000 -6.1500 
20.3100 -5.9700 
20.0000 -1.8400 
4 -40.0000 0.0000 1 
19.9900 -6.1700 
19.9900 -15.1700 
20.0100 -15.1700 
20.0100 -6.1700 
0.2 0.0000 


2.60750 0.0 0:0 
0 ‘ 
7 40.0000 0.0000 9.5000 
1.0000 20.0000 1000.0000 20.0000 13.502 
13.502 0.0000 | 

0.1000 

170.0000 170.0000 0.0000 

0.0 
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APPENDIX C 


“MATDATA” output files generated by SHIPMO.IN. These files provide regular wave 
response mass added and excitation force matrix constants for given ship, speed, and 
wave angle. File names are described in the following format: 


m -matdata. 

sork -Vessel simulated, s-SLICE; k-KAIMALINO. 

v or vh -Motion simulated, v-vertical; vh-all six degrees of freedom. 
Speed -Zero to twenty knots in one-knot increments. _ 
Angle -Zero to 180 degrees in five-degree increments. 


Example: | 
mkvh5_180.txt = Kaimalino, motion in 6-dof, at 5 knots, 180° wave angle. 


SLICE matdata files KAIMALINO matdata files 
























mkvhO_0.txt msvh0Q_0.txt 
mkvh1_0.txt msvh1- 0.txt 
mkvh2_0.txt msvh2_0.txt 
mkvh3_0.txt msvh3_0.txt 
mkvh4_0.txt msvh4_0.txt 
mkvh5_0.txt msvh5_0.txt 














msvh20_0.txt 
msvh0_5.txt 

msvh1_5.txt 
msvh2_5.txt 


mkvh20_0.txt 
mkvh0O_5.txt 
mkvh1_5.txt 
mkvh2_5.txt 









mkvh20 180.txt 


(2 ships) X (21 speeds) X (37 angles) = 1554 files 
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APPENDIX D 


% samplemain.m 
% This file takes input speed, heading, spring const., and damping constants 
% and returns regular and random wave response in the vertical plane. 
Jo Shipmo output files (matdata) are loaded. 
% Coupled heave-pitch response of the Slice and Kaimalino are modeled: 
% Parametric studies conducted on input variables. 
% Results are used in graphics program where significant force is plotted versus 
% x-coordinates specified for different values of parametric variable. 
Jo | 
% Dimensional version (U.S. units) 
% Get run info 
%o 
clear 
kcount=1; 
=35000; 
% Kaimalino resistance for V=15kts, (#) 
Ltow=[10:2:20]; 
%Range of tow lengths 
k_connection=T./Ltow; 
%HS=[0:2:10]; 
%Not used in this version since 
%HScount=1; 
%HS is defined in loop 
Vkt=15; 
Vcount=1; 
betacount=1; 
betadeg=180; 
ZYinput(‘Heading (deg) = '); 
ccount=1; 
c_connection=0; 
Yoinput;(‘Damping constant (pound.sec/ft) = '); 


while kcount~=(length(Ltow)+1), %Loop of values varied for 
% parametric study 
HScount=1; %counter for x coordinate 
HS=[0:2:10]; %X-coordinate of parametric plot — 
while HScount~=(length(HS)+1); Yocreate Y-coord for a parametric 


Zinput and range of X-coordinates. 
K_connection=k_connection(kcount)+i*c_connection(ccount); 


V_string=num2str(Vkt(Vcount)); 
beta_string=num2str(betadeg(betacount)); 
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7%The matdata output files default to the vertical only format when the 
heading angle is 0 or 180 degrees. 


ZoSet up file reading format: 
trige = 30; 
f3loc = 27; fSloc=29; 
if betadeg(betacount)== 
trigg = 27; 
f3loc = 26; f5loc=27; 
elseif betadeg(betacount)==180 
trigg = 27; 
f3loc = 26; fSloc=27; 
end 


Jo 

% Load FRONT SHIP data file msvhV_beta.txt 

Jo 

load_filename=strcat(‘msvh', V_string,'_',beta_string,'.txt’); 
filename_s=load(load_filename); 

Yo | 

% Load REAR SHIP data file 

Jo 

load_filename=strcat(‘mkvh', V_string,'_',beta_string,'.txt’); 
filename_k=load(load_filename); 


J P 
% GENERAL DATA 
Jo 
V=Vkt* 1.6878; % Convert to ft/sec 
lambda_min=20; % Min wave length (ft) 
lambda_max=1000; % Max wave length (ft) 

- delta_lambda=20; % Wave length increment (ft) 
rho=1.9905; % Water density 
zeta=1; % Regular wave height 
L=105; % Reference length 
eg=32.2; % Gravitational constant 
x_s=-46; Yo FRONT SHIP attachment point 
x_k=+40; Yo REAR SHIP attachment point 


beta=betadeg*pi/180; 
lambda=lambda_min:delta_lambda:lambda_max; 

% Vector of wavelengths 
wavenumber=2.0* pi./lambda; % Wave number 
omega=sqrt(wavenumber*g); % Wave frequency 
omegae=omega-wavenumber* V(Vcount)*cos(beta(betacount)); 

% Frequency of encounter 
period=2.0*pi./omega; 
periode=2.0*pi./omegae; 
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omega=omega’; 

omegae=omegae’; 

filesize=size(lambda); 
Jambda_size=trigg*filesize(2); 

Jo 

% FRONT SHIP 

% 

% Set mass matrix elements 

Jo 

M33s=filename_s(3:trigg:lambda_size,3); 
M35s=filename_s(3:trigg:lambda_size,5); 
M53s=filename_s(5:trigg:lambda_size,3); 
M55s=filename_s(5:trigg:lambda_size,5); 

%o 

% Added mass terms 

Jo 

A33s=filename_s(9:trigg:lambda_size,3); 
A35s=filename_s(9:trigg:lambda_size,5); 
A53s=filename_s(11:trigg:lambda_size,3); 
A55s=filename_s(11:trigg:lambda_size,5); 

Jo 

% Damping terms 

Jo 

B33s=filename_s(15:trigg:lambda_size,3); 
B35s=filename_s(15:trigg:lambda_size,5); 
B53s=filename_s(17:trigg:lambda_size,3); 
B55s=filename_s(17:trigg:lambda_size,5); 

% 

% Hydrostatic terms 

%o 

C33s=filename_s(21 :trigg:lambda_size,3); 
C35s=filename_s(21:trigg:lambda_size,5); 
C53s=filename_s(23:trigg:lambda_size,3); 
C55s=filename_s(23:trigg:lambda_size,5); 

% : 

% Total exciting forces 

Jo 
F3s_t_amp=filename_s(f3loc:trigg:lambda_size,5); 
F5s_t_amp=filename_s(f5loc:trigg:lambda_size,5); 
F3s_t_pha=filename_s(f3loc:trigg:lambda_size,6); 
F5s_t_pha=filename_s(f5loc:trigg:lambda_size,6); 
F3s_t=F3s_t_amp.*exp(i*F3s_t_pha.*pi/180.0); 
F5s_t=F5s_t_amp.*exp(i*F5s_t_pha.*pi/180.0); 
Jo 

% Froude/Krylov exciting forces 

J 


1? 








F3s_f_amp=filename_s(f3loc:trigg:lambda_size,1); 
F5s_f_amp=filename_s(f5loc:trigg:lambda_size, 1); 
F3s_f_pha=filename_s(f3loc:trigg:lambda_size,2); 
F5s_f_pha=filename_s(fSloc:trigg:lambda_size,2); 
F3s_f=F3s_f_amp.*exp(i*F3s_f_pha.*pi/180.0); 
F5s_f=F5s_f_amp.*exp(i*F5s_f_pha.*pi/180.0); 

%o 

% Diffraction exciting forces 

%o 
F3s_d_amp=filename_s(f3loc:trigg:lambda_size,3); 
F5s_d_amp=filename_s(f5loc:trigg:lambda_size,3); 
F3s_d_pha=filename_s(f3loc:trigg:lambda_size,4); 
F5s_d_pha=filename_s(f5loc:trigg:lambda_size,4); 
F3s_d=F3s_d_amp.*exp(i*F3s_d_pha.*pi/180.0); 
F3s_d=F5s_d_amp.*exp(i*F5s_d_pha.*pi/180.0); 
%o 

% REAR SHIP 

%o 

% Set mass matrix elements 

Jo 

M33k=filename_k(3:trigg:lambda_size,3); 
M35k=filename_k(3:trigg:lambda_size,5); 
M53k=filename_k(5:trigg:lambda_size,3); 
M55k=filename_k(5:trigg:lambda_size,5); 

Jo 

% Added mass terms 

% 

A33k=filename_k(9:trigg:lambda_size,3); 
A35k=filename_k(9:trigg:lambda_size,5); 
A53k=filename_k(11:trigg:lambda_size,3); 
A55k=filename_k(11:trigg:lambda_size,5); 

Jo 

% Damping terms 

%o 

B33k=filename_k(15:trigg:lambda_size,3): 
B35k=filename_k(15:trigg:lambda_size,5); 
B53k=filename_k(17:trigg:lambda_size,3); 
B55k=filename_k(17:trigg:lambda_size,5); 

Jo 

% Hydrostatic terms 

Jo 

C33k=filename_k(21:trigg:lambda_size,3); 
C35k=filename_k(21:trigg:lambda_size,5); 
C53k=filename_k(23:trigg:lambda_size,3); 
C55k=filename_k(23:trigg:lambda_size,5); 

Jo 
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% Total exciting forces 
% 
F3k_t_amp=filename_k(f3loc:trigg:lambda_size,5); 
F5k_t_amp=filename_k(f5loc:trigg:lambda_size,5); 
F3k_t_pha=filename_k(f3loc:trigg:lambda_size,6); 
F5k_t_pha=filename_k(f5loc:trigg:lambda_size,6); 
F3k_t=F3k_t_amp.*exp(i*F3k_t_pha.*pi/180.0); 
F5k_t=F5k_t_amp.*exp(i*F5k_t_pha.*pi/180.0); 
J 
% Froude/Krylov exciting forces 
% 
F3k_f_amp=filename_k(f3loc:trigg:lambda_size, 1); 
F5k_f_amp=filename_k(fSloc:trigg:lambda_size, 1); 
F3k_f_pha=filename_k(f3loc:trigg:lambda_size,2); 
F5k_f_pha=filename_k(fSloc:trigg:lambda_size,2); 
F3k_f=F3k_f_amp.*exp(i*F3k_f_pha.*pi/180.0); 
F5k_f=F5k_f_amp.*exp(i*F5k_f_pha.*pi/180.0); 
% 
% Diffraction exciting forces 

Jo | 
F3k_d_amp=filename_k(f3loc:trigg:lambda_size,3); 
F5k_d_amp=filename_k(fSloc:trigg:lambda_size,3); 
F3k_d_pha=filename_k(f3loc:trigg:lambda_size,4); 
F5k_d_pha=filename_k(f5loc:trigg:lambda_size,4); 
F3k_d=F3k_d_amp.*exp(i*F3k_d_pha.*pi/180.0); 
F5k_d=F5k_d_amp.*exp(i*F5k_d_pha.*pi/180.0); 
% | 
% MATCHING CONDITION 
Yo 
A33bar_s=-(omegae.*2).*(M33s+A33s)+i*omegae.*B33s+C33s; 
A35bar_s=-(omegae.*2).*(M35s+A35s)+i*omegae.*B35s+C35s; 
A53bar_s=-(omegae.*2).*(M53s+A53s)+i*omegae.*B53s+C53s; 
ASSbar_s=-(omegae.*2).*(M55s+A55s)+i*omegae.*B55s+C55s; 
A33bar_k=-(omegae.*2).*(M33k+A33k)+i*omegae.*B33k+C33k; 
A3Sbar_k=-(omegae.’2).*(M35k+A35k)+i*omegae.*B35k+C35k: 
A53bar_k=-(omegae.*2).*(M53k+A53k)+i*omegae.*B53k+C53k; 
A5Sbar_k=-(omegae.*2).*(M55k+A55k)+i*omegae.*B55k+C55k; 
% 
mu3_s=(A55bar_s.*F3s_t-A3Sbar_s.*F5s_t)./(A33bar_s.*A55bar_s- 
A3Sbar_s.*A53bar_s); 
nu3_s=(A55bar_s+A35bar_s*x_s)./(A33bar_s.*A55bar_s- 
A3S5bar_s.*A53bar_s); 
mu5_s=(A53bar_s.*F3s_t-A33bar_s.*F5s_t)./(A53bar_s.*A35bar_s- 
A33bar_s.*A55bar_s); 
nu5_s=(A53bar_s+A33bar_s*x_s)./(AS3bar_s.*A35bar_s- 
A33bar_s.*A55bar_s); 


79 








mu3_k=(A55bar_k.*F3k_t-A35bar_k.*F5k_t)./(A33bar_k.*A55bar_k- 
A35bar_k.*A53bar_k); 
nu3_k=(A55bar_k+A35bar_k*x_k)./(A33bar_k.*A55bar_k- 
A35bar_k.* A53bar_k); 
mu5_k=(A53bar_k.*F3k_t-A33bar_k.*F5k_t)./(A53bar_k.*A35bar_k- 
A33bar_k.*A55bar_k); 
nu5_k=(A53bar_k+A33bar_k*x_k)./(A53bar_k.*A35bar_k- 
A33bar_k.* A5Sbar_k); 

% 

a=mu3_s-mu5_s*x_s-mu3_k+mu5_k*x_k; 
b=nu3_s-nu5_s*x_s+nu3_k-nu5_k*x_k; 
f=(K_connection*a)./(1+b.*K_connection); 


J 
f_s=-f; % Connection force on FRONT SHIP 
f_k=f; % Connection force on REAR SHIP 


eta3_s=mu3_s+nu3_s.*f_s; % FRONT SHIP heave 
eta5_s=mu5_s+nu5_s.*f_s; % FRONT SHIP pitch 
eta3_k=mu3_k+nu3_k.*f_k; % REAR SHIP heave 
eta5_k=mu5_k+nu5_k.*f_k; % REAR SHIP pitch 
Xi_s=eta3_s-eta5_s*x_s; % FRONT SHIP motion at connection 
x1_k=eta3_k-eta5_k*x_k; % REAR SHIP motion at connection 
xi0_s=mu3_s-mu5_s*x_s; % FRONT SHIP motion at connection for 
% zero f 
xi0_k=mu3_k-mu5_k*x_k; % REAR SHIP motion at connection for 
% zero f 
9; | 
% Random wave calculations 
% Pierson-Moscowitz spectrum 
% 
waveheight(HScount)=HS(HScount); 
POWER =-.032*(g/HS(HScount))*2; 
S  =(0.0081*2%2).*exp(POWER./(omega.“4))./(omega.45); 
% Convert S(w) to S(we) | 
Se =S./(1-(2.0/g)*omega*V(Vcount)*cos(beta(betacount))); 
- © 
% Define response spectra 
% | 
Sf =((abs(f)).42).*Se; 
Sxi_s =((abs(xi_s)).42).*Se; 
Sxi_k =((abs(xi_k)).42).*Se; 
Sxi0_s =((abs(xi0_s)).42).*Se; 
Sxi0_k =((abs(xi0_k)).42).*Se; 
SF3s_t =((abs(F3s_t)).42).*Se; 
SF3k_t =((abs(F3k_t)).42).*Se; 
% 
% Initializations 
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% 
Sf_i=0; 
Sxi_s_1=0; 
Sx1_k_i=0; 
Sxi0_s_i=0; 
Sx10_k_i=0; 
SF3s_t_i=0; 
SF3k_t_i=0; 
Yo 
% Integral S(w)*IRAOIN2 
— GQ 
for I=2:1:filesize(2), 
Shi = Sfi + 0.5*(Sf(D + Sf(I-1)) * (omegae(I-1)- 
omegae(I)); 
Sxi_s_i = Sxi_s_i + 0.5*(Sxi_s(I) + Sxi_s(-1)) * (omegae(I-1)- 
omegae(I)); 
Sxi_k_i = Sxi_k_i + 0.5*(Sxi_k() + Sxi_k(-1)) * (omegae(I-1)- 
omegae(I)); 
Sx10_s_i= Sxi0_s_i + 0.5*(Sxi0_s(D) + Sxi0_s(J-1)) * (omegae(I- 
1)-omegae(I)); 
Sxi0_k_i= Sxi0_k_i + 0.5*(Sxi0_k(D + Sxi0_k(I-1)) * (omegae(I- 
1)-omegae(I)); 
SF3s_t_i= SF3s_t_i + 0.5*(SF3s_t() + SF3s_t(I-1)) * (omegae(I- 
1)-omegae(I)); 
SF3k_t_i= SF3k_t_i + 0.5*(SF3k_t(D + SF3k_t(1-1)) * (omegae(I- 
1)-omegae(I)); 


end 
% 
% Significant double amplitudes 
Jo 
sig f =4.0*sqrt(Sf_i); 
Sig_xi_s = 4.0*sqrt(Sxi_s_i); 
sig_xi_k = 4.0*sqrt(Sxi_k_i); 
sig_xi0_s = 4.0*sqrt(Sxi0_s_i); 
sig_x10_k = 4.0*sqrt(Sxi0_k_1); 
sig_F3s_t = 4.0*sqrt(SF3s_t_i); 
sig _F3k_t = 4.0*sqrt(SF3k_t_i); 
sig_fk(HScount,kcount)=sig_f; %o(xcoord-row, param. var-col) 
HScount=HScount+1; % X-coord. counter 
end 
kcount=kcount+1; %Parametric variable counter 
end 
% 
ZYocall graphics program: 
samplegraphpm 


$samplegraphpm.m 
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%This program is for output graphics from samplemain.m 

%oParametric plots of significant force vs. x-variable plotted 

%ofor varied parametric values. 

format bank 

figure(1) %Sig. force Amplitudes vs HS for 
Yovaried k's (T/Ltow) 

kstrng=num2str(k_connection); 

Jowavestr=num2str(waveheight); 

cstmg=num7?2str(c onan), 

plot(HS, sig_tk) 

grid 

titlstr=[' Sig. force vs. H_{1/3} for \beta= ',beta_string,'‘o',’ V= ',V_string,'kts ',' 

c=',cstrng]; | 

title([titlstr]) 

xlabel("H_{ 1/3 }') 

ylabel(\sigma_f’) 
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Jomaingraph.m 

%This program is for output graphics from samplemain.m 

% 

figure(1) 

%absolute motion magnitude vs omega 

kstrng=num2str(k_connection); 

cstrng=num2str(c_connection); 
plot(omega,abs(xi_s),'S',omega,abs(xi_k),'o',omega,abs(xi_s),omega,abs(xi_k)) 
grid 

titlstr=['Absolute motion magnitude vs frequency (rad/sec) for k= ‘,kstrng,'’ c=’, 
cstmg]; 

title((titIstr]) 

xlabel(‘omega’) 

ylabelCMAGNITUDE(x1_s,xi_k)') 

stringl='xi_s'; string2='xi_k’; 

legend(string1,string2) 

%J ; 

% 

figure(2) 

Yabsolute motion phase angle vs omega 
plot(omega,57.32*angle(xi_s),'S',omega,57.32*angle(xi_k),'o'’,omega,57.32*angle 
(x1_s),omega,57.32*angle(xi_k)) 

grid | 
titlstr2=['Absolute motion phase angle vs frequency (rad/sec) for k= ',kstrng,' c=", 
cstmg]; 

title([titlstr2]) 

xlabel(‘omega’) 

ylabel(PHASE ANGLE(xi_s,xi_k)’) 

legend(string1,string2) 

% 

% | 

figure(3) 

%Connection force magnitude vs omega 
plot(omega,abs(f_s),'S',omega,abs(f_s)) 

grid | 
titlstr=[' Connection force magnitude vs frequency (rad/sec) for k=',kstrng, 
c=', cstrng]; | 

title([titlstr]) 

xlabel(‘omega’) 

ylabel(‘Connection force magnitude’) 

string1='f_s’; 

legend(string1) 

%o 


% 
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figure(4) 

Connection force phase angle vs omega 
plot(omega,57.32*angle(f_s),'S',omega,57.32*angle(f_s)) 
grid 

titlstr=['Connection force phase angle vs frequency (rad/sec) for k=',kstrng,’ c=", 
cstrng]; 

title([titlstr]) 

xlabel(‘omega’) 

ylabel(‘Connection force phase angle’) 

string 1='f_s'; 

legend(string1) 
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